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List  of  Figures 

1  A  schematic  depicting  the  vertically -resolved  one-dimensional  coupled  canopy  model  (solid  black 

lines).  The  vertical  distribution  of  canopy  density  and  associated  percentage  of  each  grid  volume 
occupied  by  sunlit  versus  shaded  leaves  and  their  scattering/absorbing  efficiency  for  photosynthet- 
ically  active  and  near-infrared  radiation  is  prescribed.  In  addition,  soil  moisture  and  temperature 
profiles  and  incoming  solar  radiafion  are  imposed.  Af  each  afmospheric  grid  volume  resolving 
fhe  vegefafion,  a  leaf-energy  balance  is  solved  for  sunlif  and  shaded  leaves  based  upon  fhe  fofal 
radiafion  absorbed  (assuming  a  random  disfribufion  of  nine  differenf  leaf  angles  for  scattering  and 
absorbing  direcf  and  diffuse  radiafion)  and  fhe  local  afmospheric  winds,  femperafure,  and  mois- 
fure.  The  exchange  of  wafer  vapor  and  heaf  befween  fhe  leaves  and  fhe  surrounding  afmosphere 
is  Iherefore  confrolled  fhrough  a  balance  of  afmospheric  demand,  phofosynfhesis,  and  soil  wafer 
availabilify  fo  fhe  frees  roofs  fhrough  boundary  layer  and  sfomafal  conducfances  for  moisfure  and 
heal.  In  fhe  furbulence-resolving  large  eddy  simulafion  code,  a  separafe  implemenfafion  of  fhis 
ID  canopy-resolving  land-surface  model  resides  af  each  horizonfal  grid  poinf  and  exfends  verfi- 
cally  from  fhe  surface  info  fhe  afmosphere  up  lo  canopy  lop.  Currenfly,  fen  grid  poinls  resolve  fhe 
lOm-fall  canopy  and  four  grid  poinls  resolve  fhe  Im-deep  soil .  6 

2  Verlical  profiles  of  fhe  mean  wind  normalized  by  fhe  friction  velocily  af  canopy  lop  {u) / n*  (lefl), 
polenfial  femperafure  minus  fhe  canopy-lop  polenlial  femperafure  in  Kelvin  (middle-lefl),  verlical 
velocily  variance  normalized  by  u\  (middle-righl),  and  vertical  velocily  skewness  (righf)  during 
CHATS  when  Ihere  were  leaves  on  fhe  frees.  The  profiles  are  crealed  by  averaging  over  fhe  same 
lime-periods  during  which  fhe  dispersion  experimenls  look  place,  which  are  1-6  PDT,  12-17  PDT, 

and  18-24  PDT  on  June  9,  2007.  The  dashed  line  represenls  fhe  lop  of  fhe  canopy,  z//i  =  1.  ...  7 

3  Horizonlally-  and  lime-averaged  verlical  profiles  of  fhe  mean  wind  normalized  by  fhe  friclion 

velocily  af  canopy  fop  (lefl),  polenfial  femperafure  minus  fhe  canopy-lop  polenlial  lem- 

peralure  in  Kelvin  (middle-lefl),  vertical  velocily  variance  normalized  by  (middle-righl),  and 
vertical  velocily  skewness  (righf)  from  fhe  coupled-canopy  furbulence-resolving  large  eddy  sim¬ 
ulations  mimicking  fhe  periods  during  CHATS  when  fhe  dispersion  experimenls  look  place.  The 
profiles  are  calculafed  from  Iwo  nearly  idenlical  simulalions  lhal  only  vary  in  fhe  geoslrophic 
wind  driving  fhe  flow;  Ug  =  (8,  4)  ms^^  for  cases  (LES-WS,  LES-S),  respectively,  where  Ihese 
geoslrophic  wind  variafions  resulf  in  Iwo  differenf  slabilily  condilions:  weakly  slable  (EES-WS) 


and  slable  (EES-S).  Nole  lhal  fhe  axis-ranges  differ  slighlly  befween  Eigures  2  and  3 .  7 

4  Horizonlally-  and  lime-averaged  vertical  profiles  of  fhe  percenfage  of  fhe  heal  flux  (w'9')  faking 

place  in  each  of  fhe  four  quadranls  (Q1-Q4,  see  insef)  for  fhe  Iwo  EES  cases;  EES-WS  (lefl)  and 
EES-S  (righf).  Nole  lhal  for  any  variable  %+  signifies  x'  >  0,  and  %  signifies  x'  <  0 .  8 

5  Verlical  profiles  of  SEg  concenlralion  af  x  =  40  m  downslream  of  a  line  source  showing  fhe  diurnal 
variafion  of  fhe  profiles  during  CHATS;  dala  are  from  June  9,  2007  wilh  hour  of  observation  noled 

nexl  lo  lines .  10 
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1.  Introduction 


Dispersion  in  the  stable  boundary  layer  (SBL)  is  important  for  a  range  of  problems  ineluding  air  quality  and 
toxie  gas  releases,  but  existing  information  on  SBL  dispersion  is  defieient  due  to  our  limited  understanding  of 
SBL  turbulenee.  This  defieieney  is  attributed  to  the  weakness  of  the  turbulenee  and  the  range  of  proeesses 
eomplieating  the  SBL  physies.  Large-eddy  simulation  (LES)  has  proven  to  be  a  powerful  tool  in  simulating 
turbulenee  in  the  planetary  boundary  layer  (PBL)  not  only  for  the  oft-studied  eonveetive  PBL  but  more  reeently 
for  the  stable  PBL.  Thus  far,  the  SBL  studies  have  been  aimed  mostly  at  weakly  stable  boundary  layers  over 
horizontally  homogeneous  surfaees.  In  addition,  the  LES  veloeity  fields  have  been  used  to  simulate  dispersion  in 
these  PBE  flows  by  traeking  “passive”  partieles  using  a  Eagrangian  partiele  dispersion  model  (EPDM). 

Under  previous  support  from  the  Army  Researeh  Offiee  (ARO),  we  developed  EPDMs  for  dispersion  in 
eonveetive  and  stable  PBEs  ineluding  some  with  a  forest  eanopy  and  showed  that  predieted  eoneentration  fields 
agreed  well  with  laboratory  data  and  field  observations  (e.g.,  Weil  et  al,  2004,  2006,  2012).  Eor  the  ARO  program 
just  eompleted,  we  extended  the  previous  effort  to  investigate  dispersion  in  more  stable  PBEs  and  SBEs  over  hor¬ 
izontally  heterogeneous  surfaees.  This  effort  was  pursued  in  part  using  our  new  eoupled  multi-layer  eanopy  -  soil 
model  {Patton  et  al,  2008).  The  program  eonsisted  of  three  main  investigations.  The  first  (Seetion  2)  had  two 
parts:  a)  further  development  of  the  eoupled-eanopy  EES  and  eomparison  with  observations  from  the  CHATS 
(Canopy  Horizontal  Array  Turbulenee  Study)  field  program  {Patton  et  al,  201 1),  and  b)  modeling  of  eanopy  dis¬ 
persion  using  the  EPDM-EES  approaeh  and  assessment  of  the  modeled  dispersion  against  the  observed  CHATS 
dispersion  data.  The  seeond  (Seetion  3)  investigation  targeted  the  effeets  of  surfaee  temperature  heterogeneity  on 
dispersion  in  the  SBE.  The  third  study  a)  used  the  eoupled  model  to  produee  a  more  stable  boundary  layer  and 
applied  the  EES  fields  generated  from  this  to  drive  our  EPDM  (Seetion  2),  and  b)  developed  an  EPDM  based  on 
parameterized  turbulenee  profiles  for  modeling  dispersion  in  more  stable  PBEs  (Seetion  4). 

2.  Coupled-Canopy  LES  of  the  PBL  and  Dispersion  Modeling 

2.1  Coupled-Canopy  LES  and  CHATS  Meteorological  Data  Analysis 

Under  earlier  ARO  funding,  we  developed  a  eoupled  eanopy-soil  model  by  extending  the  NOAH  land-surfaee 
system  to  inelude  a  vertieal  distribution  of  the  sensible-,  radiative-,  and  latent-heat  fluxes,  a  height-dependent 
leaf-area  density,  stomatal  resistanee  and  leaf  boundary  layer  resistanees.  This  new  land-surfaee  model  is  suit¬ 
able  for  eoupling  with  either  one-dimensional  models  of  eanopy  turbulenee  and  sealar  transport  that  eould  be 
implemented  in  mesoseale  models,  or  with  turbulenee  resolving  simulations  like  those  generated  by  large-eddy 
simulation  (Pig.  1).  The  intent  is  to  allow  for  dynamie  ealeulation  of  vegetation-imposed  sealar  sources/sinks 
determined  by  atmospherie  demand,  and  to  provide  an  avenue  for  the  fluid  meehanies  to  respond  to  rapidly- 
fluetuating  leaf  temperatures.  Under  our  eurrent  funding,  this  eanopy-resolving  land-surfaee  model  has  been 
eoupled  to  NCAR’s  large-eddy  simulation  eode.  This  eoupled  modeling  system  serves  as  the  basis  for  generating 
the  eanopy-influeneed  wind  fields  in  our  CHATS-foeused  investigations  of  stability  influenees  on  contaminant 

Table  1:  Bulk  parameters  characterizing  the  two  LES  cases,  LES-WS:  Ug  =  8  ms^\  and  LES-S:  Ug  =  4  ms^K 
The  quantities  presented  are:  1)  the  friction  velocity  at  canopy  top  h,  where  m*  =  {{u'w')^  +  2)  the  PBL 

depth,  Zi  defined  as  the  height  of  the  maximum  vertical  gradient  in  virtual  potential  temperature  above  the  canopy,  3) 
the  Monin-Obukhov  length  evaluated  at  canopy  top,  L  =  — [g/(0y)]  (w'0y)),  and  4)  the  bulk  Richardson  number 
Rib  =  {[g/  (0v)]  A0y  Az)/([A(m)]^  -|-  [A(v)]^)  evaluated  over  two  different  height  ranges  Az;  from0<  z<h  and  from  h<z<  Zt- 


Uif  Zi  L  h/L  Zi/L  Rib  Rig 

(ms^^)  (m)  (m)  0<z<h  h<z<Zi 

LES-WS  0.29  172  204  0.05  0.85  0.02  0.13 

LES-S  0.08  106  12  0.85  9.01  1.39  0.31 
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Figure  1:  A  schematic  depicting  the  verticaiiy-resoived  one-dimensionai  coupied  canopy  modei  (soiid  biack  iines).  The 
verticai  distribution  of  canopy  density  and  associated  percentage  of  each  grid  voiume  occupied  by  suniit  versus  shaded 
ieaves  and  their  scattering/absorbing  efficiency  for  photosyntheticaiiy  active  and  near-infrared  radiation  is  prescribed.  In 
addition,  soii  moisture  and  temperature  profiies  and  incoming  soiar  radiation  are  imposed.  At  each  atmospheric  grid  voiume 
resoiving  the  vegetation,  a  ieaf-energy  baiance  is  soived  for  suniit  and  shaded  ieaves  based  upon  the  totai  radiation  absorbed 
(assuming  a  random  distribution  of  nine  different  ieaf  angies  for  scattering  and  absorbing  direct  and  diffuse  radiation)  and 
the  iocai  atmospheric  winds,  temperature,  and  moisture.  The  exchange  of  water  vapor  and  heat  between  the  ieaves  and  the 
surrounding  atmosphere  is  therefore  controiied  through  a  baiance  of  atmospheric  demand,  photosynthesis,  and  soii  water 
avaiiabiiity  to  the  trees  roots  through  boundary  iayer  and  stomatai  conductances  for  moisture  and  heat.  In  the  turbuience- 
resoiving  iarge  eddy  simuiation  code,  a  separate  impiementation  of  this  ID  canopy-resoiving  iand-surface  modei  resides  at 
each  horizontai  grid  point  and  extends  verticaiiy  from  the  surface  into  the  atmosphere  up  to  canopy  top.  Currentiy  ten  grid 
points  resoive  the  1  Om-taii  canopy  and  four  grid  points  resoive  the  Im-deep  soii. 


dispersion. 

Under  Task  la,  simulations  using  the  coupled  LES-canopy  system  have  been  used  to  investigate  stably  strat¬ 
ified  canopy  turbulence.  These  simulations  used  512  x  512  x  256  grid  points  to  resolve  a  512  x  512  x  256  m^ 
domain.  The  simulations  were  generated  by  initiating  the  calculations  with  weak  solar  forcing  and  allowing  the 
sun  to  set.  As  the  solar  forcing  diminishes,  the  coupled  canopy  slowly  cools  and  subsequently  cools  the  atmo¬ 
sphere  in  contact  with  it.  Two  different  stabilities  were  generated  by  modifying  the  geostrophic  wind  (Ug),  where 
Ug  =  4  or  8  ms^^  Table  1  presents  bulk  parameters  more  fully  characterizing  the  simulations. 

Task  la  also  focused  on  evaluating  the  coupled  canopy  model  across  a  wide  range  of  atmospheric  stability  and 
establishing  the  mechanisms  responsible  for  turbulent  exchange  between  the  canopy  layers  and  those  aloft  with 
increasing  importance  of  buoyancy  forces.  Towards  that  end,  data  analysis  of  the  CHATS  30-m  meteorological 
tower  was  performed  (Dupont  and  Patton,  2012a,b)  to  investigate  the  variation  of  velocity,  temperature,  and 
humidity  fields  and  their  statistical  properties  with  atmospheric  stability  and  seasonal  canopy  changes,  and  the 
mechanisms  leading  to  momentum/scalar  transport.  This  CHATS  data  analysis  serves  to  evaluate  the  simulations. 

The  CHATS  meteorological  high-rate  data  from  the  30-m  tower  revealed  the  sensitivity  of  velocity,  tempera¬ 
ture  and  humidity  fields  within  and  above  the  walnut  canopy  to  atmospheric  stability  and  seasonal  canopy  changes 
(Dupont  and  Patton,  2012a,b).  This  analysis  was  performed  both  prior  to  and  following  leaf-out  of  the  trees  and 
hence  resulted  in  both  de-foliated  and  foliated  periods.  Profiles  of  meteorological  statistical  properties  were  an¬ 
alyzed  for  five  stability  regimes  and  two  seasonal  periods.  In  addition,  an  analysis  in  terms  of  flow  structures 
within  and  above  the  canopy  was  made  to  highlight  the  flow  character  using  the  well-established  mixing-layer 
analogy  (Finnigan,  2000;  Finnigan  et  al,  2009).  In  near-neutral  conditions,  the  canopy  flow  better  resembles  a 
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Figure  2:  Vertical  profiles  of  the  mean  wind  normalized  by  the  friction  velocity  at  canopy  top  (u)  / m*  (left),  potential 
temperature  minus  the  canopy-top  potential  temperature  in  Kelvin  (middle-left),  vertical  velocity  variance  normalized  by  ul 
(middle-right),  and  vertical  velocity  skewness  (right)  during  CHATS  when  there  were  leaves  on  the  trees.  The  profiles  are 
created  by  averaging  over  the  same  time-periods  during  which  the  dispersion  experiments  took  place,  which  are  1-6  PDT, 
12-17  PDT,  and  18-24  PDT  on  June  9,  2007.  The  dashed  line  represents  the  top  of  the  canopy,  z/h  =  \. 


Figure  3:  Horizontally-  and  time-averaged  vertical  profiles  of  the  mean  wind  normalized  by  the  friction  velocity  at  canopy 
top  (m) /m*  (left),  potential  temperature  minus  the  canopy-top  potential  temperature  in  Kelvin  (middle-left),  vertical  veloc¬ 
ity  variance  normalized  by  (middle-right),  and  vertical  velocity  skewness  (right)  from  the  coupled-canopy  turbulence¬ 
resolving  large  eddy  simulations  mimicking  the  periods  during  CHATS  when  the  dispersion  experiments  took  place.  The 
profiles  are  calculated  from  two  nearly  identical  simulations  that  only  vary  in  the  geostrophic  wind  driving  the  flow;  Ug  =  (8, 
4)  ms^^  for  cases  (LES-WS,  LES-S),  respectively,  where  these  geostrophic  wind  variations  result  in  two  different  stability 
conditions:  weakly  stable  (LES-WS)  and  stable  (LES-S).  Note  that  the  axis-ranges  differ  slightly  between  Eigures  2  and  3. 
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Figure  4:  Horizontally-  and  time-averaged  vertical  profiles 
of  the  percentage  of  the  heat  flux  (w'9')  taking  place  in  each 
of  the  four  quadrants  (Q1-Q4,  see  inset)  for  the  two  LES 
cases;  LES-WS  (left)  and  LESS  (right).  Note  that  for  any 
variable  %:  signifies  >  0,  and  signifies  x!  <  0. 


plane  mixing-layer  during  the  foliated  period  than  in  the  defoliated  period,  whereas  in  the  no-leaf  period  the  flow 
appears  more  as  the  superposition  of  a  wall  boundary-layer  (e.g.,over  a  flat  plate)  and  a  plane  mixing-layer  flow. 
The  plane  mixing-layer  analogy  weakens  with  deviations  from  near-neutral  stability,  whether  unstable  or  stable. 
For  example,  in  unstable  conditions  the  importance  of  mixing-layer  type  structures  in  the  canopy  decreases  as 
exhibited  by  weaker  canopy-induced  shear  of  streamwise  velocity,  lower  velocity  skewness  and  kurtosis  maxima, 
and  reduced  efficiency  of  the  turbulent  transport  of  momentum.  In  addition,  with  greater  instability,  heat  and 
momentum  appear  to  be  transported  by  canopy-scale  thermal  plumes  instead  of  the  mixing-layer  type  structures. 

Figure  2  presents  profiles  of  the  mean  wind  normalized  by  the  canopy-top  friction  velocity  (m)/m*,  poten¬ 
tial  temperature  deviation  (9)  —  (9)/,  (in  K),  dimensionless  vertical  velocity  variance  and  vertical  velocity 

skewness  from  CHATS  averaged  over  the  five-six  hour  time  periods  on  June  9  during  which  the  disper¬ 

sion  experiments  took  place  (1-6,  12-17,  18-24)  PDT.  The  simulated  mean  wind  profiles  under  stable  conditions 
(Fig.  3)  exhibit  shapes  similar  to  those  observed  at  CHATS  (Fig.  2).  Like  those  at  CHATS,  the  stably-stratified 
simulations  reveal  significant  variation  across  stability  regimes,  where  the  weak  wind  case  develops  strong  strat¬ 
ification  across  the  canopy  top  (middle  panel.  Fig.  3)  and  the  horizontal  winds  diminish  in  the  subcanopy  due  to 
decoupling  with  the  overlying  flow.  The  potential  temperature  profile  for  the  strongly  stable  case  (middle  panel) 
also  mimics  the  CHATS  profiles  for  similar  stability;  the  temperature  inversion  near  the  canopy  top  is  important 
in  limiting  the  vertical  dispersion  of  scalars  emitted  within  the  canopy.  Another  significant  feature  is  the  sign 
change  of  the  vertical  velocity  skewness  with  increasing  stability;  a  feature  also  found  in  the  CHATS  data  (Fig. 
2).  This  positive  within-canopy  vertical  velocity  skewness  indicates  that  vertical  motions  in  the  subcanopy  are 
dominated  by  small-scale  within-canopy  buoyant  plumes  rather  than  by  canopy-induced  mixing-layer  eddies  and 
has  significant  implications  for  canopy  dispersion. 

Quadrant  analysis  of  scalar  fluxes  decomposes  the  flux  into  quadrants  based  upon  the  sign  of  the  fluctuating 
quantities  contributing  to  the  co-variance  (e.g.,  Willmarth  and  Lu,  1972)  and  helps  quantify  the  character  of  the 
turbulent  motions  performing  the  transport.  Fig.  4  presents  a  quadrant  analysis  based  upon  horizontal-  and  time- 
averages  of  the  two  stably-stratified  large  eddy  simulations.  It  is  important  to  note  that  Fig.  4  presents  the  quadrant 
analysis  in  magnitude  fraction  form  (i.e.,  the  absolute  value  of  the  flux  in  a  particular  quadrant  normalized  by  the 
sum  of  the  absolute  value  of  the  flux  across  all  four  quadrants)  which  eliminates  the  sign  of  the  flux  and  forces 
the  sum  over  all  four  quadrants  to  a  value  of  one. 

For  weakly  stable  conditions  (LES-WS,  left  panel  Fig.  4),  the  majority  of  the  flux  occurs  in  quadrants  Q2 
(downward  moving  warm  air)  and  Q4  (upward  moving  cool  air),  with  Q2  becoming  increasingly  important  in 
the  canopy’s  vicinity  (i.e.,  ,  below  z/h  =  2).  This  heat  flux  partitioning  is  consistent  with  organized  turbulent 
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motions  generated  by  veloeity  shear  at  the  eanopy-top  transporting  warm  air  downward  into  the  eanopy  from  aloft 
penetrating  nearly  to  the  ground  surfaee  with  slightly  weaker  importanee  of  eool  within-eanopy  air  transported 
upwards  out  of  the  eanopy  (e.g.,  Chen,  1990;  Katul  et  al,  1997;  Thomas  and  Token,  2007).  With  inereasing 
stability  (LES-S,  right  panel  Fig.  4),  Q2  and  Q4  events  still  dominate  transport  in  the  upper  eanopy  and  above, 
but  do  not  reveal  as  signifieant  an  inerease  in  Q2  events  with  deseent  toward  the  eanopy  top  eompared  to  LES- 
WS.  Deep  in  the  eanopy  {z/h  <  ~0.6),  relatively  warm  soils  beneath  the  relatively  eool  eanopy-top  air  generates 
vertieal  motions  sueh  that  these  layers  are  dominated  by  upward  moving  warm  air  (Ql)  and  downward  moving 
eool  air  (Q3).  The  leaves  in  the  eanopy’s  upper  reaehes  are  exposed  to  the  sky  in  both  simulations,  however 
in  EES-S,  these  leaves  radiatively  eool  more  quiekly  than  the  rate  at  whieh  warm  air  is  transported  downward 
from  aloft  ereating  the  relatively  strong  stratifieation  oeeurring  at  z//i  =  ~0.8  (Fig.  3,  seeond  panel).  This 
strongly  stratihed  layer  aets  like  a  barrier  to  vertieal  transport,  thereby  suggesting  that  the  shear-generated  eddies 
at  eanopy-top  do  not  (or  infrequently)  penetrate  all  the  way  to  the  surfaee;  a  feature  signifieantly  impaeting  the 
eseape  effieieney  and/or  vertieal  exehange  of  eontaminants  emitted  within  the  eanopy’s  lower  levels. 

2.2  CHATS  Dispersion  Data  Analysis  and  LPDM-LES  Modeling 

The  main  purpose  of  the  CHATS  dispersion  experiment  was  to  determine  the  effeets  of  eanopy-indueed  stability 
on  vertieal  dispersion,  espeeially  at  night  when  air  in  the  lower  eanopy  was  unstable  but  air  in  the  upper  eanopy 
and  above  it  was  stable.  This  was  intended  to  fill  a  void  in  eanopy  dispersion  experiments,  i.e.,  most  eanopy 
dispersion  experiments  have  been  eondueted  during  near-neutral  eonditions.  For  the  6-day  experiment,  there 
were  two  daytime  periods  (of  about  4-6  hours  duration)  and  five  nighffime  periods  of  a  similar  durafion.  In 
fhe  experimenf,  an  SFg  fraeer  was  released  from  a  100-m  long  line  souree  in  fhe  10-m  deep  walnuf  orehard;  fhe 
souree  was  loeafed  40  m  upwind  of  fhe  30-m  fall  mefeorologieal  lower  and  was  elevaled  aboul  a  meler  above  fhe 
surfaee.  The  key  ambienl  measuremenls  were  SFg  eoneenlralions  oblained  from  seven  heighls  on  fhe  lower  and 
nine  ground-level  SFg  samplers  loeafed  af  dislanees  of  20  m  lo  120  m  from  fhe  souree. 

Slabilily  effeels  on  dispersion  over  fhe  diurnal  eyele  were  well-demonslraled  by  verlieal  profiles  of  fhe  hourly- 
averaged  SFg  from  fhe  30-m  lower.  During  daylime  (Fig.  5a),  fhe  observalions  showed  lhal  fhe  eoneenlralion 
deereased  upwards,  oflen  wilh  a  eoneave  upwards  shape,  whieh  is  lypieal  of  a  surfaee  release  in  an  unslable  or 
eonveelive  boundary  layer.  In  Fig.  5a,  fhe  above-eanopy  slabilily  ranged  from  weakly  unslable  (early  afternoon) 
lo  weakly  slable  in  Ihe  late  afternoon,  while  exhibiling  weak  (slabilily,  inslabilily)  in  Ihe  (upper,  lower)  eanopy 
layers,  respeelively.  However,  al  nighl  Ihe  profiles  exhibited  mueh  higher  surfaee  eoneenlralions  and  a  more  rapid 
deerease  wilh  heighl  al  leasl  for  Ihe  higher  wind  eases  before  midnighl  (Fig.  5b).  The  kink  in  Ihe  SFg  profile 
(z  =  7  m;  Fig.  5b)  near  Ihe  plume  lop  is  in  middle  of  an  elevated  inversion  or  slable  layer  in  Ihe  upper  pari  of 
Ihe  eanopy;  e.g., see  Ihe  polenlial  lemperalure  profile  for  Ihe  slable  FES  ease  (FES-S)  in  Fig.  3.  Due  lo  Ihe  slable 
slralifiealion,  verlieal  mixing  would  be  suppressed  al  shorl  Iravel  limes  or  dislanees  (i.e.,  v  =  40  m.  Fig.  5)  bul 
exhibil  a  more  gradual  heighl  varialion  furlher  downslream  as  shown  in  Fig.  8. 

In  Ihe  weak  and  variable  winds  after  midnighl  and  in  Ihe  early  morning  (17/,  ~  0.4  ms^^),  Ihe  surfaee  eon¬ 
eenlralions  were  lypieally  higher  and  more  variable,  and  Ihe  verlieal  SFg  profiles  tended  towards  fuller  shapes, 
somelimes  eoneave  downwards  (Fig.  5e).  In  Ihis  Ihird  period,  Ihe  slabilily  as  represented  by  Ihe  verlieal  polenlial 
lemperalure  gradienl  was  nol  as  large  as  during  Ihe  pre-midnighl  period,  whieh  may  explain  Ihe  somewhal  greater 
plume  deplh  in  Fig.  5e  lhan  in  Fig.  5b. 

In  our  FPDM  approaeh,  we  Iraek  single  parlieles  to  oblain  Ihe  mean  dispersion  and  eoneenlralion  field.  When 
using  FES  veloeilies  to  drive  Ihe  model,  we  divide  Ihe  parliele  veloeily  into  “resolved”  and  “sublilter-seale”  (SFS) 
eomponenls,  eonsislenl  wilh  Ihe  FES.  The  SFS  veloeily  is  a  random  value  oblained  from  a  sloehaslie  model  using 
Ihe  SFS  lurbulenl  kinelie  energy.  In  addilion,  Ihe  FPDM-FES  approaeh  ean  be  used  to  oblain  an  ensemble  of 
dispersion  realizalions  for  determining  Ihe  mean,  rool-mean-square  devialion,  and  fluelualing  fields  of  dispersion 
quanlilies;  Ihis  has  been  done  for  Ihe  eonveelive  boundary  layer  (Weil  et  al,  2012). 

The  FES  fields  for  Ihe  weakly-slable  (FES-WS)  and  slable  (FES-S)  simulations  diseussed  in  Seetion  2.1  were 
used  to  obtain  the  FPDM  mean  eoneentration  fields;  aboul  1.1  x  10^  parlieles  were  released  from  Ihe  modeled 
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SFg  Concentration  (ppt) 


Figure  5:  Vertical  profiles  of  SF^  concentration  atx  =  40  m  downstream  of  a  line  source  showing  the  diurnal  variation  of 
the  profiles  during  CHATS;  data  are  from  June  9,  2007  with  hour  of  observation  noted  next  to  lines. 


line  source  in  the  lower  canopy.  The  LPDM  vertical  concentration  profile  is  compared  with  the  CHATS  SFg  data 
for  the  daytime  cases  (Fig.  6)  and  has  a  shape  generally  similar  to  the  measurements,  but  underestimates  the 
near-surface  concentrations  by  about  20%  on  average.  This  is  likely  due  to  the  stability  differences — modeled 
versus  observed — over  the  6-hr  period.  However,  other  comparisons  (Figs.  7  and  8)  with  the  surface  sampler  data 
shows  quite  good  agreement,  and  thus  there  is  the  possibility  of  small  differences  between  the  ground-level  and 
tower  measurements. 

The  surface  concentration  and  its  variation  with  downwind  distance  is  a  key  quantity  of  interest.  LPDM 
calculations  were  made  for  SBL  conditions  that  approximated  but  were  not  identical  to  the  field  observations, 
fhus  requiring  scaling  of  fhe  resulfs.  For  similar  mefeorological  condifions,  fhe  concenfrafion  scaling  for  a  canopy 
line  source  of  sfrengfh  Qi  is  CUhh/Q(,  =  /(2f/,),  where  X/,  =  [Uhh)  is  a  dimensionless  time  or  disfance,  x  is 
fhe  disfance  from  fhe  source,  and  C  is  fhe  concenfrafion.  The  dimensionless  concenfrafion  versus  X/,  fhen  should 
be  fhe  same  in  fhe  model  and  field  and  was  used  fo  scale  fhe  LPDM  runs  for  fhe  LES-WS  and  LES-S  fields  fo 
fwo  observafional  periods:  1)  fhe  near-neufral  daytime  case  (Eig.  5a,  which  for  fhe  sake  of  comparison  will  be 
called  WS),  and  2)  fhe  early  nighttime  case  (Eig.  5b,  and  will  be  called  S).  The  wind  and  furbulence  fields  from 
EES-WS  are  good  approximations  for  WS,  and  fhose  from  EES-S  are  only  rough  approximafions  for  S. 

Eigure  7  shows  fhe  concenfrafion  versus  disfance  for  fhe  fwo  periods.  Eor  fhe  WS  case  (green  line),  fhe  EPDM 
resulf  agrees  well  wifh  fhe  observations,  and  fhe  concenfrafion  falls  off  approximafely  as  C  oc  near  fhe  source 
(x  <  40  m),  which  is  consisfenf  wifh  dispersion  in  a  PBE  wifh  no  canopy  {e.g.,  van  Ulden,  1978).  However,  for 
40  m  <  X  <  120  m,  fhe  concenfrafion  decreases  more  sharply  as  C  oc  x^^,  which  has  no  counferparf  for  fhe  PBE 
wifhouf  a  canopy.  This  decrease  resulfs  from  bofh  furbulenf  dispersion  and  fhe  canopy-induced  sfrong  verfical 
shear  of  fhe  horizonfal  wind  speed  across  fhe  canopy  fop  (leff  panel,  Eigs.  2  and  3).  Over  fhe  40  m  fo  120  m 
disfance,  fhe  mean  plume  height  Zp  increases  from  O.Th  to  l.5h  with  a  distance  dependence  like  Zp  x,  and  the 
vertical  plume  spread  exhibits  a  similar  dependence  (not  shown).  Eurthermore,  for  this  height  change,  the 
mean  wind  U  increases  by  a  factor  of  three  (Eig.  3),  and  the  mean  plume  speed  varies  as  Up  oc  x.  Since  the 
concentration  behaves  as  C  oc  Qi/ {UpG^),  Up{x)  provides  for  additional  dilution  of  the  plume  as  shown  in  Eig.  7, 
i.e.,  to  produce  the  C  oc  x^^  dependence.  Eor  greater  distances  {i.e.,  x  >  120  m),  the  concentration  tends  towards 
a  C  oc  x^^/^  variation,  which  is  consistent  with  the  long-time  concentration  behavior  in  a  PBE  with  no  canopy. 

The  EPDM  results  for  the  stable  (S)  case,  scaled  to  the  conditions  for  the  nighttime  period  in  Pig.  5b,  are 
given  by  the  blue  lines  in  Pig.  7.  Due  to  the  very  weak  canopy-top  winds  iUh  =  0.48  ms^'  in  the  EES),  the 
EPDM  mean  concentration  field  is  developing  wifh  fime,  and  resulfs  are  shown  for  fwo  emission  fime  durations, 
38  min  and  64  min,  wifh  fhe  laffer  exhibifing  high  concenfrafions  exfending  furfher  downwind.  This  frend  of 
increasing  concenfrafion  wifh  disfance  and  fime  is  expecfed  fo  confinue,  fhus  requiring  longer-fime  simulations. 
Neverfheless,  fhe  S  case  (64  min)  exhibifs  higher  concenfrafions  fhan  fhe  WS  case  as  would  be  expected,  and  fhe 
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Figure  6:  LPDM  and  observed  vertical  profiles  of  Figure  7;  LPDM  and  observed  SFe  surface  concentration  as  a 
SF6  concentration  downstream  of  a  line  source  during  function  of  distance  downwind  of  a  line  source  during  CHATS 
daytime  in  CHATS;  data  from  June  9,  2007.  on  June  9,  2007.  Points  and  error  bars  show  the  mean  and  root- 

mean-square  deviation  in  concentration  over  each  measurement 
period. 


model  results  agree  approximately  with  the  observations  at  very  short  range,  x  <  40  m.  We  believe  that  model 
predictions  for  longer  emission  times  would  approximate  better  the  measurements  at  greater  distances. 

Figure  8  shows  the  modeled  time  or  distance  evolution  of  the  vertical  profiles  of  concentration  for  the  two 
LPDM-LES  cases  and  the  important  differences  between  them.  The  LPDM-LES  results  from  the  EES-WS  case 
exhibit  a  “lapse”  or  continuously  decreasing  vertical  gradient  of  the  dimensionless  concentration  with  z/h  and  Xh 
as  a  result  of  the  weakly-stable  stratification  within  and  above  the  canopy.  In  contrast,  EPDM-EES  results  using 
data  from  the  EES-S  case  leads  to  an  elevated  within-canopy  scalar  concentration  inversion  with  approximately 
uniform  (or  slightly  increasing)  values  below  the  inversion;  the  uniform  value  results  from  the  weak  convection 
and  turbulence  below  the  inversion  (Eig.  3)  which  generate  vertical  mixing  driven  by  the  contrast  between  the 
relatively  warm  soil  and  the  leaf-induced  cool  air  in  the  upper  canopy.  Essentially,  the  plume  mimics  the  thermal 
structure  but  the  scalar  inversion  strength  continually  decreases  with  X/,  due  to  the:  1)  small  “line”  nature  and 
dimension  of  the  scalar  source  versus  the  larger  areal  nature  of  the  thermal  source,  and  2)  the  stronger  wind  and 
turbulence  above  the  canopy-induced  inversion  which  “mixes  out”  the  scalar  from  the  canopy  and  decreases  its 
inversion  strength  with  increasing  time  or  X/,. 

3.  Surface  Temperature  Heterogeneity  Effects  on  Dispersion  in  the  SBL 

This  study  follows  previous  work  by  Stoll  and  Porte -Agel  (2009)  and  Mironov  and  Sullivan  (2010)  who  showed 
how  the  surface-temperature  heterogeneity  modified  the  SEE  structure  and  turbulent  transport  properties  relative 
to  a  surface-temperature  homogeneous  SEE.  Eoth  boundary  layers  were  cases  of  a  weakly  stable  boundary  layer 
(WSEE)  with  moderate-to-strong  winds.  In  Mironov  and  Sullivan  (2010),  the  temperature  heterogeneity  was 
generated  with  spanwise  homogeneous  surface  temperature  strips  that  had  a  sinusoidal  variation  in  the  x  (mean 
wind  («))  direction  and  a  horizontal-mean  surface  temperature  that  was  the  same  as  in  the  homogeneous  SEE. 
Eoth  studies  showed  that  the  heterogeneous  (HET)  SEE  was  more  turbulent  and  well-mixed  with  respect  to 
temperature  than  the  homogeneous  (HOM)  SEE.  An  outstanding  finding  from  Mironov  and  Sullivan  (2010)  was 
that  the  temperature  variance  peaked  both  at  the  surface  and  at  the  SEE-top  with  maximum  values  far  greater 
than  in  the  HOM  case.  Throughout  the  bulk  of  the  HET  SEE,  the  temperature  variance  was  relatively  uniform  or 
“well-mixed”  which  contrasted  with  the  more  vertically-inhomogeneous  structure  of  the  HOM  SEE. 


11 


CU^h/Q, 


Figure  8:  Evolution  of  LPDM  vertical  profiles  of  the  dimensionless  concentration  with  dimensionless  distance  X/,  down¬ 
stream  of  a  line  source  in  weakly-stable  (WS)  and  stable  (S)  boundary  layers.  LPDM  WS  results  compared  with  CHATS 
surface-sampler  data  for  daytime  (hours  12-  17  on  June  9,  2007);  points  and  error  bars  have  definition  as  in  Fig.  7. 


We  used  the  LES  fields  from  Mironov  and  Sullivan  (2010)  to  drive  the  LPDM  and  investigate  the  dispersion 
characteristics  for  a  continuous  source  at  different  heights  in  the  two  SBLs.  Heterogeneity  effects  can  be  assessed 
by  comparison  to  our  earlier  results  for  the  homogeneous  SBL  (Weil  et  al,  2006).  This  comparison  was  simplified 
since  the  LES  fields  in  the  Mironov  and  Sullivan  (2010)  work  were  created  for  the  same  domain  size  [(400  m)^], 
number  of  grid  points  (200  x  200  x  192),  geostrophic  wind  Ug  (=  8  ms^^),  initial  stratification,  and  other  condi¬ 
tions  as  for  the  HOM  SBL  described  in  Weil  et  al.  (2006).  Lor  the  HEX  SBL,  the  Zi  =  221  m,  n*  =  0.27  ms^', 
and  the  vertically-averaged  wind  speed  U  was  7.4  ms^^  which  are  all  close  to  their  counterparts  for  the  HOM 
case. 

We  conducted  LPDM  simulations  for  several  source  heights,  but  our  main  focus  was  on  a  surface  release  in 
order  to  explain  the  large  surface  temperature  variance  found  by  Mironov  and  Sullivan  (2010).  Calculations  were 
made  for  sources  at  the  leading  edge  of  the  hot  and  cold  strips  and  in  the  middle  of  each  strip,  where  the  strip 
length  in  the  x  direction  was  200  m;  hence,  sources  were  located  at  =  0,  100,  200,  and  300  m  along  the  x 
axis  with  about  17  sources  spanning  the  strip  in  the  y  direction.  Lor  each  y-strip  of  sources,  10  realizations  of 
the  concentration  field  were  obtained  using  a  sequence  of  initialization  times  that  were  about  1  min  apart,  and 
1.28  X  10^  particles  were  released  in  each  realization. 

Our  results  showed  an  oscillatory  behavior  of  the  mean  concentrations  with  distance  downwind  of  the  source. 
This  was  attributed  to  the  sinusoidally- varying  surface  temperature  which  generated  gravity  waves  in  this  stratified 
flow.  An  understanding  of  this  behavior  can  be  made  with  reference  to  the  Sullivan  and  McWilliams  (2002)  study 
of  turbulent  flow  over  water  waves  in  the  presence  of  stable  stratification,  where  the  surface  waves  generated 
gravity  waves  in  the  air  above  the  water.  Lollowing  their  analysis,  we  decomposed  the  instantaneous  concentration 
c  as  c{x,y,z,t)  =  C{x,y,z)  -\-Cw{x,z)  -|- c'(x,y,z,f),  where  C  is  the  ensemble-mean  concentration,  is  a  phase- 
averaged  concentration,  c'  is  a  turbulent  fluctuation,  and  x,  y,  z  are  measured  with  respect  to  the  source  location. 
The  ensemble-mean  concentration  for  a  particular  line  or  strip  of  sources  (along  y  at  a  fixed  x^)  is  (c)  =  C  -|-  c^, 
where  (  )  denotes  an  ensemble  average.  Our  focus  is  on  C  and  c„. 

The  mean  concentration  C  was  obtained  by  averaging  the  concentrations  from  a  large  number  (289)  of 
sources  evenly  distributed  over  the  LES  x-y  domain  as  in  Weil  et  al.  (2004).  Here,  we  analyze  the  crosswind- 
integrated  concentration  (CWIC)  O'  nondimensionalized  as  C^Uzi/ 2  as  a  function  of  the  dimensionless  distance 
X  =  M*x/ {Uzi),  where  Q  is  the  source  strength.  The  LPDM  surface  CWIC  results  for  the  HET  and  HOM  SBLs 
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Figure  9:  Dimensionless  crosswind-integrated  integrated  concentration  (CWIC)  at  the  surface  as  a  function  of  dimension¬ 
less  distance  due  to  a  surface  source  in  a  surface-temperature  heterogeneous  SBL:  a)  source  (x^  =0)  at  leading  edge  of  hot 
strip,  b)  source  (xs  =  Lx/l)  at  leading  edge  of  cold  strip,  and  c)  source  fc  =  LxjA-)  in  middle  of  hot  strip.  denotes  either 
the  mean  CWIC  O'  or  the  average  CWIC  over  a  y-strip  of  sources  at  fixed  Xs,  i.e.,  (c^).  Ensemble-mean  CWIC  from  the 
heterogeneous  and  homogeneous  SBLs  shown  for  reference.  Here,  Xj  =  and  Xj  =  Lx. 


were  in  excellent  agreement  with  one  another  as  shown  in  Fig.  9.  This  assures  one  of  the  similarity  in  the  vertical 
structure  of  the  mean  wind  and  turbulent  velocity  statistics  in  these  two  SBLs;  these  variables  drive  the  mean 
dispersion.  The  mean  CWIC  results  also  agree  well  with  field  observations  (not  shown)  from  the  Prairie  Grass 
experiments  {Barad,  1958). 

Figure  9a  also  shows  the  mean  dimensionless  CWIC  (c^)  versus  X  for  sources  at  the  leading  edge  of  the 
hot  strip  (Xi  =  0).  The  (c'’)  is  the  average  of  the  10  realizations,  which  do  not  exhibit  a  large  scatter,  and  only 
the  realizations  are  shown;  the  actual  average  will  be  presented  in  the  future.  As  is  clear,  the  (c^)  oscillates 
over  the  entire  distance  range  shown  which  is  about  three  times  the  domain  length  Lx  =  400  m  or  =  0.064 
in  dimensionless  distance  (X)  units.  Note  that  the  (c^)  oscillates  with  a  wavelength  X„  that  is  about  half  of  the 
surface  temperature  wavelength  Xj,  i.e.,  Xw  =  Xj /2.  The  shorter  wavelength  for  the  c„  can  be  explained  by 
reference  to  the  Sullivan  and  McWilliams  (2002)  study,  which  showed  that  the  wave-induced  heat  flux  O^w^ 
generated  by  the  water  waves  had  a  wavelength  one-half  of  the  water  surface  wavelength.  This  occurred  because 
the  wave-induced  temperature  (On,)  and  vertical  velocity  (w^)  perturbations  were  90°  out-of-phase.  We  believe 
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that  the  same  phenomena  exists  in  this  heterogeneous  situation  and  explains  the  shorter  Xw 

Figure  9b  shows  similar  results  (red  lines)  for  sourees  loeated  along  the  leading  edge  of  the  hot  strip  {xs  = 
Lxjl  =  200  m)  and  overlaid  on  those  (blue  lines)  from  Fig.  9a.  The  averages  or  (c^)  versus  X  from  these  two 
sourees  are  in-phase,  whieh  further  supports  the  shorter  wavelength  Xw  for  the  €„  field.  Finally,  Fig.  9e  presents 
results  for  a  souree  at  the  midpoint  of  the  hot  strip  {xs  =  Xy^jA),  and  shows  that  the  (c^)  is  90°  out-of-phase  with 
the  results  for  x^  =  0. 

The  net  result  for  from  a  distribution  of  surfaee  sourees  over  the  x  —  y  domain  of  a  HET  SBL  is  that  the 
ensemble-mean  O  eoneentration  would  be  similar  to  the  HOM  SBL  as  shown  above,  but  the  varianee  would 
not.  The  varianee  would  eontain  a  c^(x,z)  term  applieable  at  distanee  x  downwind  of  a  partieular  souree,  and 
for  a  distribution  of  sourees,  these  varianee  terms  would  be  superposed  to  produee  a  total  varianee  greater  than 
the  that  due  to  the  SBL  turbulenee  alone,  i.e.,  {c'^).  In  the  HET  SBL  analyzed  by  Mironov  and  Sullivan  (2010), 
the  temperature  varianee  was  quite  large  implying  a  non-negligible  eontribution  from  or  due  to  the  surfaee 
temperature  distribution. 

4.  A  Parameterized  LPDM  for  Moderate  to  Very  Stable  Boundary  Layers 

This  study  aimed  at  the  development  of  a  relatively  simple  model  for  understanding  and  predieting  the  behavior  of 
the  mean  wind,  potential  temperature,  and  flux  profile  behavior  in  fhe  SBL  over  a  broad  range  of  sfabilify — from  a 
weakly-sfable  (WSBL)  fo  a  very-sfable  (VSBL)  boundary  layer.  The  model  was  based  on  fhe  equations  governing 
fhe  above  variables  and  used  an  eddy-diffusivily  (K)  elosure  from  fhe  Brost  and  Wyngaard  (1978)  model.  This 
was  infended  fo  provide  some  preliminary  informalion  for  moderately-  fo  very-sfable  PBLs,  where  LES  eurrenfly 
does  nof.  The  modeled  or  parameferized  profiles  of  fhe  mean  wind,  femperafure,  and  furbulenee  for  fhese  more 
sfable  eases  ean  be  and  have  been  used  fo  drive  an  LPDM  for  dispersion  prediefions  in  fhese  eondifions.  We  label 
Ibis  fhe  LPDM-PT,  whieh  is  based  on  parameferized  furbulenee  (PT)  and  winds. 

An  eddy-diffusivily  {K)  approaeh  is  juslified  for  modeling  fhe  mean  wind,  femperafure,  ele.  in  fhe  SBL 
beeause  of  fhe  small  furbulenee  lengfh  seales  in  fhe  SBL  as  found  from  field  dafa  and  LES.  In  addition  fo  use 
of  fhe  Brost  and  Wyngaard  (1978)  K  profiles,  fhere  are  fwo  ofher  imporfanf  aspeefs  of  fhis  ‘W-model”  approaeh 
{Weil,  2011):  a)  reduefion  of  fhe  model  equafions  fo  a  simple  ordinary  differenlial  form  using  a  similarify  variable, 
and  b)  modeling  of  fhe  flux  profiles  using  a  simple  K  approaeh.  The  model  has  been  fesfed  in  fwo  SBL  limifs 
and  agrees  well  wifh:  1)  LES  resulfs  and  field  dafa  for  a  WSBL  {zi/L  =  1.6),  and  2)  field  observations  for  a 
VSBL  {Zi/L  =  42),  where  Zi  and  L  are  fhe  SBL  heighf  and  Monin-Obukhov  lengfh,  respeelively.  Note  lhal  in  fhe 
following  figures  and  diseussion,  we  also  use  h  =  Zi  to  denote  the  SBL  height. 

Ligure  10  presents  the  evolution  of  the  modeled  temperature  differenee  profiles  for  a  VSBL,  z,/L  =  42,  and 
eompares  fhem  wifh  fhe  CASES99  field  dafa  {Mahrt  and  Vickers,  2006).  The  agreemenf  is  good  and  shows  fhaf 
fhe  SBL  dimensionless  profiles  do  nof  ehange  mueh  wifh  lime.  This  is  due  fo  fhe  very  small  diffusivily  in  fhe 
VSBL  by  eomparison  fo  lhal  in  fhe  WSBL,  where  fhe  larger  K  leads  fo  a  greater  lime  varialion  and  differenl 
profile  shapes. 

Ligure  1 1  presenfs  fhe  ‘W-model”  SBL  wind  profiles  for  a  broad  range  of  sfabilify,  where  fhe  sfabilify  index 
Z,/L=  1.6,  10,  and  40.  The  modeled  profiles  agree  approximately  wifh  ofher  resulfs  in  fwo  eomparisons:  fhe  LES 
predielion  (dashed  blue  line)  for  fhe  WSBL  {Zi/L  =  1.6)  and  fhe  CASES99  field  dafa  for  a  VSBL  {Zi/L  =  42), 
whieh  oeeurred  for  a  weak  wind  Uh  (—  3  m/s)  al  fhe  SBL  lop.  An  allernalive  sfabilify  index  is  fhe  bulk  Riehardson 
number  Ribu  whieh  is  based  on  fhe  SBL  deplh,  fhe  polenlial  femperafure  ehange  A0/,  over  fhe  SBL,  and  fhe  Uh- 
Lor  CASES99,  fhe  Ribh  —  1.2  and  is  al  or  near  fhe  erifieal  value  or  limil  {i.e.,  1)  for  furbulenee  fo  exisl  or  be 
mainlained.  A  eomparison  for  mueh  weaker  winds  {Uh  —  0.35  m/s)  in  fhe  LLOSS  experimenl  {Mahrt  and  Vickers, 
2006)  showed  lhal  fhe  model  did  nof  work  well  for  fhe  wind  and  flux  profiles;  here,  fhe  Ribh  was  51  and  henee 
mueh  greater  lhan  1.  Af  LLOSS,  fhe  Uh  was  so  low  lhal  meehanieally-foreed  furbulenee  was  exlremely  weak  and 
dominafed  by  ofher  foreings,  e.g.,densily-driven  flows  and  furbulenee.  As  a  resull,  we  lenfafively  believe  lhal  fhe 
“operating  range”  of  fhe  LPDM  based  on  fhe  parameferized  winds  and  furbulenee  is  for  Zi/L  <  40  or  perhaps 
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Figure  10:  Profiles  of  dimensionless  potential  temperature 
as  a  function  of  time  for  a  VSBL,  ZilL  =  42;  here  h  =  Zi,  the 
SBL  height. 


Figure  11:  Mean  wind  profile  in  the  stable  boundary  layer  from  the 
“K-model”  as  a  function  of  the  stability  index  h/L  for  hjh  =1.6 
(green  line),  10  (red  line),  and  40  (blue  line).  K-model  prediction 
can  be  compared  with  LES  result  for  a  WSBL  (h/L=  1.6)  and 
CASES99  data  for  a  VSBL  (h /L  =  42). 


somewhat  larger  and/or <  1.2.  These  limits  need  further  clarification. 

For  concentration  estimates  using  the  LPDM-PT,  we  could  adopt  the  A'-model  wind  profile  buf  have  inifially 
chosen  a  simpler  mefhod  in  which  fhe  MO  profile  is  correcfed  for  fhe  sfabilify.  The  correcfion  is  an  exponenfially- 
decaying  funcfion  of  z/zt,  exp(— az/z;),  where  fhe  paramefer  a  is  empirically  chosen  and  depends  on  Zi/L.  Figure 
12  shows  fhaf  fhe  correcfed  MO  profile  is  a  good  mafch  fo  fhe  .K'-model  profile  over  mosf  of  fhe  SBL.  In  addifion, 
fhe  LPDM-PT  requires  fhe  vertical  furbulence  velocify  Ow  Figure  13  shows  a  parameferizafion  (red  line) 
compared  wifh  a  range  of  resulfs:  a)  Nieuwstadfs  (1984)  field  dafa,  LES  resulfs,  and  wind-funnel  dafa  all  for 
a  WSBL  wifh  Zi/L  ~  1.6,  and  b)  fhe  CASES99  field  dafa  for  Zi/L  ~  42.  We  believe  fhaf  fhe  parameferizafion 
is  an  adequafe  fif  fo  fhese  diverse  resulfs  alfhough  further  work  and  dafa  are  necessary  fo  fesf  or  improve  fhe 
parameferizafion  for  more  sfable  condifions,  i.e.,  fhe  VSBE. 

Predicfions  of  fhe  crosswind-infegrafed  concenfrafion  (CWIC,  C^)  due  fo  a  surface  source  in  fhe  SBE  have 
been  made  using  fhe  EPDM-PT  for  a  range  of  sfabilify.  The  model  was  compared  wifh  surface  CWIC  observafions 
from  fhe  Prairie  Grass  experimenf  (Barad,  1958),  which  was  conducted  over  a  grass  surface  and  for  disfances  x  < 
800  m.  Eigures  14  and  15  show  fhe  dimensionless  CWIC,  C^Uzi/Q,  as  a  function  of  fhe  dimensionless  disfance 
X  from  bofh  fhe  model  and  experimenfs  for  fwo  groups  of  sfabilify,  which  range  from  weakly-fo-moderafely 
sfable.  Alfhough  fhere  is  a  some  bias  and  scaffer,  fhe  frend  and  overall  agreemenf  befween  fhe  model  and  dafa  is 
considered  quife  good.  Eurfher  examination  of  fhe  Prairie  Grass  dafa  record  and  ofher  experimenfs  will  be  made 
in  an  affempf  fo  obfain  measuremenfs  for  more  sfable  condifions,  Zi/L  >  5,  for  evaluating  fhe  model. 

5.  Summary 

This  research  program  extended  our  previous  sfudies  of  furbulence  and  dispersion  in  fhe  planefary  boundary 
layer  (PBE)  fo  more  sfable  PBEs  including  fhose  inferacfing  wifh  a  foresl  canopy,  and  fo  sfable  boundary  layers 
(SBEs)  over  horizonfally-heferogeneous  surfaces.  There  were  four  key  findings.  Eirsf,  our  coupled-canopy  large- 
eddy  simulafion  (EES)  model  was  a)  pushed  fo  more  sfable  PBEs  fhan  previously  achieved  wifh  simulation 
of  fhe  vertical  sfrucfure  of  mean  and  furbulence  quanfifies  wifhin  and  above  fhe  canopy,  including  one  case 
producing  a  1  K  femperafure  inversion  in  fhe  upper  canopy,  and  b)  verificalion  of  fhis  temperature  structure  by 
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U/u. 

Figure  12:  Comparison  between  the  “corrected”  Monin- 
Obukhov  (MO)  wind  profiles  and  the  K-model  profiles  in 
the  SBL  for  weakly-  to  moderately-stable  conditions;  here, 
M*  is  the  surface  friction  velocity 


Figure  13:  A  parameterization  of  the  dimensionless  ver¬ 
tical  turbulence  component,  Owju^,,  and  comparison  with 
field  data,  wind-tunnel  data,  and  LES  results  for  a  WSBL 
(hlL=  1.6)  and  CASES99  field  observations  fora  VSBL. 


comparison  with  field  observations  in  weakly-  to  strongly-stable  conditions  in  an  orchard  canopy.  Second,  the 
LES  fields  were  used  to  drive  our  Lagrangian  particle  dispersion  model  (LPDM)  for  concentration  predictions  that 
a)  matched  observations  in  the  above  two  stabiliy  regimes,  b)  showed  perhaps  for  the  first  time  the  importance 
of  the  strong  canopy-top  wind  shear  in  accelerating  the  decrease  in  surface  concentrations  with  downstream 
distance,  and  c)  produced  vertical  profiles  of  scalar  concentration  that  ranged  from  “lapse”  conditions  to  one 
with  a  scalar  inversion  for  strong  stability.  Third,  the  LPDM  investigations  of  dispersion  in  an  SBL  with  either  a 
heterogeneous  or  homogeneous  surface-temperature  distribution  showed  important  effects  of  the  heterogeneity. 
Although  the  mean  concentrations  in  the  HET  and  HOM  cases  were  similar,  heterogeneity  can  produce  signficant 
surface  concentration  variance  resulting  from  heterogeneity-induced  gravity  wave-induced  scalar  flux;  this  finding 
has  the  potential  for  determining  or  identifying  the  degree  of  surface-temperature  heterogeneity  in  the  SBL  from 
in-situ  measurements,  satellite  observations,  and  modeling.  Lourth,  modeled  or  parameterized  profiles  of  wind, 
temperature,  and  turbulence  were  derived  for  weakly-  to  very-stable  PBLs  and  provide  the  inputs  for  driving 
an  LPDM  based  on  parameterized  turbulence  (PT)  for  stronger  stratification.  The  LPDM-PT  was  found  to  give 
predictions  of  surface  concentrations  that  matched  observations  in  weakly-  to  moderately-stable  PBLs. 
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observations  (points). 
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